Tensor to scalar ratio of perturbation amplitudes and inflaton dynamics 
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For the inflaton perturbations it is shown that the evolution of the difference between the spectral 
indices can be translated into information on the scale dependence of the tensor to scalar amplitudes 
ratio, r, and how the scalar field potential can be derived from that information. Examples are given 
where r converges to a constant value during inflation but dynamics are rather different from the 
power-law model. Cases are presented where a constant r is not characteristic of the inflationary 
dynamics though the resulting perturbation spectra are consistent with the CMB and LSS data. 

The inflaton potential corresponding to r given by a n-th order polynomial of the e-folds number 
is derived in quadratures expressions. Since the observable difference between the spectral indices 
evaluated at a pivot scale yields information about the linear term of that polynomial, the first 
order case is explicitly written down. The solutions show features beyond the exponential form 
corresponding to power-law inflation and can be matched with current observational data. 

PACS numbers: PACS numbers: 98.80.Cq, 98.80.Es, 98.70.Vc 



I. INTRODUCTION 

During the early phase of the cosmological evolution 
the universe could experience a period of accelerated ex- 
pansion known as inflation. The simplest scenario with 
a time dependent equation of state yielding the required 
negative pressure is that of a single real scalar field, the 
inflaton, with dynamics dominated by its potential en- 
ergy. For extensive details and references on this scenario 
and on the other topics mentioned in this Introduction 
see book § and Ref. @. 

The inflaton quantum fluctuations would be stretched 
by the expansion beyond the radius within which causal 
interactions take place. These curvature perturbations 
could reenter the causal horizon in much later epochs 
and, through gravitational collapse, lead to anisotropies 
in the cosmic microwave background (CMB) tempera- 
ture and to cosmological large scale structure (LSS). Dur- 
ing inflation primordial gravitational waves are also pro- 
duced. These tensor perturbations induce a curled po- 
larization in the CMB radiation and increase the overall 
amplitude of its anisotropies at large angular scales. 

A hint about the physics in the very early universe 
could be obtained by fitting to data the result of analyt- 
ical calculations of the CMB and density spectra. These 
calculations depend on the values of some parameters 
which are the ones to be pinned down. The initial con- 
ditions for the evolution of the thermal anisotropies are 
also characterized by several quantities. These primor- 
dial parameters are often given as the multipoles of 
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where A stands for the normalized amplitudes of the 
scalar (As) or tensor (At) perturbations, the corre- 
sponding spectral indices, n, are defined by, 
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k = aH is the comoving wavenumber corresponding to 
the wavelength matching the Hubble horizon during in- 
flation and £;* is a pivot scale. 

Initially, only As(k*) was fitted. Since inflation pre- 
dicts nearly scale-invariant spectra, the tilt given by the 
scalar spectral index was then taken into account. This 
reduces the primordial spectrum to a power-law function 
of the scale k. The only single field model exactly yielding 
such a spectrum is power-law inflation jj| where the cos- 
mic scale factor, a, behaves like a power-law of the cosmic 
time, t, and the inflaton potential is an exponential func- 
tion. Signals of nonzero 'running' of the scalar index, 
dn$ /dink, have already been reported Q and must be 
refined by near future observations, allowing, this way, 
to move beyond the power-law approximation. 

The role of the tensor perturbations deserves also at- 
tention when determining the best-fit values of the cos- 
mological parameters from CMB and LSS spectra. That 
is motivated in part by the possibility of measuring the 
cosmic background polarization, allowing the tensorial 
contribution to be indirectly determined. This contribu- 
tion can be parametrized in terms of the relative ampli- 
tudes of the tensor and scalar perturbations, 
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where a is a constant. Presently, due to measurement 
limitations, a constant value of r is fitted. 

Nevertheless, in Ref. it was shown that few infla- 
tionary models produce an exactly constant tensor to 
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scalar ratio, and, in order to be proper scenarios of in- 
flation, they must be observationally indistinguishable 
from power-law inflation, where r = constant too. Since 
power-law inflation is just one of many suitable final 
stages of the inflaton dynamics ra], it implies that, ei- 
ther a constant value for r may be no characteristic of 
the underlying inflationary evolution or that conclusions 
about the inflaton potential beyond an exponential form 
may be no possible to be drawn |Q. This last factor 
could be a strong limitation for programs of the inflaton 
potential reconstruction. 

Taking the above into account, the aim of this letter 
is to analyze the relation between the functional form of 
the tensor to scalar ratio and the inflationary dynamics. 
After introducing in Sec. ID the relevant equations, the 
analysis is done in Sec. pT[| by means of several exam- 
ples. Finally, the results are discussed in Sec. IV and 
it is concluded that the difference between the spectral 
indices can be a very useful quantity because it yields in- 
formation on the scale dependence of the tensor to scalar 
ratio, hence allowing to observe features of the inflaton 
potential different from the exponential form. 



II. KEY EQUATIONS 

The horizon-flow functions proved to be useful quanti- 
ties while describing the inflationary dynamics and per- 
turbations. They were defined in Ref. Jq] as 



MN) 
dm 



£m+l 



d In | e„ 
dN 



m > , (5) 



where djj = a/ a denotes the Hubble distance, N = 
ln(a/ai) the number of e-folds since time t{, and dm = 
dn(ti). During inflation, < t\ < 1. For m > 1, e m 
may take any real value, though usually |e m | < 1 for the 
inflationary spectra to be weakly scale dependent. For 
instance, in Ref. || the CMB constrains t\ < 0.1 and 
|e 2 | < 0.3 were found for fixed e 3 = 0. 

To next-to-leading order in an expansion in terms of 
the horizon-flow functions the spectral indices (||) and 
(0) can be written as || [u| 

n s -l = -2e 1 -e 2 -2e?-(2C + 3)e 1 e 2 -Ce 2 e3,(6) 
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where C « —0.7293. If these expressions are analyzed as 
a system of differential equations (where the differentia- 
tion, according with definitions (||) , is done with respect 
to N), then it is not a closed system. Even after differ- 
entiating Eq. (R) and substituting the result and Eq. (0) 
itself into Eq. (0) , one still has to deal with three indepen- 
dent variables, namely, ej., ns and ut (see Refs. || [ll[] 
for more detailed discussion). There arc not other inde- 
pendent equations relating these variables which could be 
used to close the system of equations (^|) and (]?]) . There- 
fore, information on the functional forms of the observ- 
ables ns and hit it is necessary in order to describe the 



dynamics of ei. This is the basic philosophy behind the 
Stewart-Lyth inverse problem |llf|. 

It follows from definitions (@)T(§) and (@) that 

dhur . , „. 

= An = n T - (n s - 1) . (8) 

This way, any information on the evolution of both spec- 
tral indices can be used as information on the scale de- 
pendence of the tensor to scalar ratio. Using expression 
(|), after subtracting Eq. (§) to Eq. (g) yields 
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Then, changing variables in Eq. (JsJ) it is obtained, 
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where Eqs. (||) were used and terms with order higher 
than the second one were neglected consistently with the 
approximations applied to derive Eqs. (^) and (g) JTo[] . 

Taking into account the definition (||) for e 2 , Eq. (|11|) 
can be rewritten as, 
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where ro is an integration constant. Now, the number of 
variables to solve for has been reduced to two by com- 
prising in r the information on ns and tit- Remains only 
one equation to deal with thus, it is still compulsory to 
find out how r behaves to describe the dynamics of e%. 
According with definitions (0) , 
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where H is an integration constant. On the other hand 
using the definition (|5|) for t\ and the Friedmann equation 
for the inflaton cosmology pi fL2|, 
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dN 

where <p is the inflaton, n 
Planck mass. Hence, given a solution for £i(N) the cor- 
responding potential as function of the inflaton field is, 



87r/mpj and mpi is the 
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V(N) = §■ [3 - ei] exp ( -2 / e^N ) , 



where the expression for the potential as function of N 
is derived from the Einstein equations for the scalar field 
cosmology |l], 11 . In this way, the inflaton potential 
can, in principle, be determined from just the information 
on the functional form of the tensor to scalar ratio or, 
equivalently, on the evolution of the difference between 
the tensor and scalar spectral indices. 
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III. SOME SOLUTIONS 

Current analysis of CMB and LSS observations yields 
r < 1 0, hence solutions satisfying lnr < shall be 
searched for. Eq. Jl2] ) can be written as, 



r ei exp(Ce 2 ) 
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Expanding this expression to next-to-leading order and 
comparing with the corresponding expansion for At /As 
|Tof it can be determined that r = 16 [|[ |||. For large 
values of |ea|, models with e 2 > are favored. For the 
next-to-leading order approximation to proceed it must 
be required the minimum value, e 2 = —0.3. Together 
with the constrain r < 1, this lead to t\ < 0.05. These 
values coincides with the 2a bounds reported in Ref . M . 

Both, Eqs. (11) and (|l2| ) are in the class of first order 
linear differential equations like, 



y' + ay = g{x) . 



(17) 



with a prime denoting differentiation with respect to x 
and a a constant. Using the ansatz y = exp(— ax)f(x), 
the general solution of Eq. ( p7[ ) is obtained as, 



y = exp(— ax) 
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where B is the constant arising from the integration of 
the homogeneous equation and the integration limits are 
determined by the range of x where g{x) is continuous. 
Thus, one has the option of solving any of Eqs. ([□]) or 
(|l2|), regarding the complexity of the corresponding left 
hand side expression. If g(x) — > /3 = const, while x — > oo, 
then, for a > 0, applying the L'Hopital's rule it can 
be determined that, for any lower integration limit, the 
solution converges to 13/ a. However, for a < the only 
solution converging to (3 /a is, 



y = — exp(— ax) / g(x) exp(ax)dx . 
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bmce m Eqs. (||) and (H), a = 1/C < 0, this implies 
an attractor given by the particular solution t\ = r*/16 
for r = = const, (or An = from Eq. (B)). This 
corresponds to power-law inflation and it is just a special 
case (that for B = 0), even if r — ► r* while N — ► oo. For 
instance, if it is assumed 



lnr = 6cxp(— aN) + lnr c 
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with Too, a > and b some constants then, according 
with Eq. (pj[), the asymptotic solution for ei will be 



ei = eioo exp 



Ca - 1 



exp(— aN) 
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where e^ = r^/16. This solution converges to t\ = 
e loo while N — > oo. For b > (b < 0), this happen 
starting from values larger (smaller) than eioo- If b = 0, 



corresponding to r = = const., then the outcome 
is just the power-law solution given by ei(AT) = eioo- 
This is the kind of behavior observed in Ref. (l3) for 
models with constant tensorial index and also in the first 
case analyzed in Ref. (l2) for models with weakly scale 
dependent spectral indices. 

According to Eq. (19), the more general solution with 
r given by Eq. (EOT) is, 
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Now the dynamics is different. If B > (B < 0) the 
solution increases (decreases) very fast starting from, 



(23) 



ei = eioo exp (B) exp - 



Ca-1 



If b = 0, this behavior is modified only in the rate of 
change of ei , consistently with the solution in Ref. [|| for 
models with r — const.. Thus, B is a bifurcation param- 
eter for the inflationary dynamics described by Eq. (p^). 

Functional forms of r converging to a constant value 
are motivated by difficulties of measuring the primordial 
gravitational waves background fjl|. Maybe the unique 
for a while alternative for determining the contribution 
of tensor modes to CMB is the derivation of the tensor to 
scalar ratio from the background radiation polarization, 
the more optimistic expectations being measuring a con- 
stant central value of r 0. However, there are other be- 
haviors of r consistent with current constrains on the pri- 
mordial spectra. For example, if in Eq. (|2(]) a is changed 
for —a, the general solution (22) is modified in exactly 
the same way, i.e., a for —a. Nevertheless, now we have 
richer dynamics as for r as for ei, depending on the val- 
ues of the different parameters involved. For B = and 
b < 0, r — * with N — > oo but there are two different 
behaviors for ei; if < a < —1/C, t\ — > oo, else ei — > 0. 
Interesting is the case B = and b > 0, since r — > oo 
with N — > oo but for < a < — 1/C, ei — > 0. According 
to Eq. (^|) , this example can be made consistent with cur- 
rent analysis of CMB and LSS observations ^ ^, |9| by 
tunning the values of a and b to control the behavior of 
e 2 and €3. For B ^ more possibilities arise depending 
on whether B is positive or negative, and greater or less 
than 6/(1 + aC). Solutions with r — * 00, ei — > and 
Nl < 1, I e 3 1 < 1 while N — > 00, can be easily found. 

Equally satisfactory, from the observational point of 
view, are functional forms of r exhibiting a transient pe- 
riod of almost power-law behavior during a sufficiently 
large numbers of N. For example, a proper choice of the 
constants a, b, a, (3 ensures that to be the case for 

r = r {at&n[a(N - N )]-b} 

(3 



x exp 



x{N - N )} (atan[a(2V - N )} - b) 



and the corresponding solution for ei, 

ei = £i(o) {atan[a(iV - N )} - b} 



(24) 
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This result is similar to case 2 in Ref. for models with 
weakly scale dependent spectral indices. 

Now, if g(x) in Eq. (|l7|) is given precise enough by, 
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with n some finite integer, then solution (|T^) reads, 
y(x) — Bq exp(— ax) 



p=0 q=0 v ^ H > 
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for p> q. The solution of Eq. ( |l2| ) is then 
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where 6 



p. 9 



a p (-1) 9 C(« +1 )p !/(p - q) !. The general 
expression for the inflaton potential corresponding to r 
given by a n-th order polynomial of the e-folds number, 
N, can be derived in terms of quadratures expressions 
for ei, by substituting solution ( p6[ ) into Eq. jl5|). 

Unfortunately, the future prospects for obtaining in- 
formation from observations for higher order coefficients 
in expansion ( |25| ) are strongly limited. However, it is 
timely to recall that dr/dhxk = rAn which, according to 
Eqs. (||) and ([!(]) is a next -to-leading order expression, 
hence a\ « aoAn(No). With these regards, it becomes 
important to analyze in details case, 



In — = a + cn{N - N ) , 
16 

where the corresponding solution for t\ is, 



ei(N) = e 1(i) exp 



Bexp 



AT 
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exp ( A/V) 
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(28) 



with ex(j) = £i(o) exp{[ao + (A^o — C) ai] C} and A = 
aiC. The asymptotes of this solution for B ^ will 
be mainly determined by the value and sign of B. How- 
ever, in the same fashion as it was shown in Ref. [g for 
A = (i.e., lnr/16 = a ), if the model yielding ei given 
by Eq. fl2g ) is expected to be compatible with current 
data, B has to be chosen an extraordinarily small num- 
ber, making the corresponding scenario very difficult to 
be observationally distinguished from power-law infla- 
tion. More interesting is the case B = 0. From Eq. (15), 
the corresponding inflaton potential is obtained as, 
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where ip = n/2{4> + tfi ). 
tential are plotted for A = 
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FIG. 1: Sectors of the inflaton potential given by Eq. (^9|) for 
A = -0.0073 (Vr) and A = 0.0073 (Vu). 




FIG. 2: The spectral indices for Vi and Vu with eiu) — 0.05 



(Vr/). In Fig. || are presented the corresponding to this 
model scalar and tensorial spectral indices, 



n s = 1 - A + [2 + (2C + 3)A] £ 



ir 



-ei(i) 



2£ 



-ei(i) 



(30) 



n T = 2 [1 + (C + 1)A] £ w 

A' 



2£ 2 



(31) 



where Eqs. @, (0) and® were used and £v^[x] is the 
Lambert TF function |1J]. For a large set of A and eui) 
values, these spectra agree with CMB and LSS observa- 
tions, n s = 1.033 ± 0.66, n T = 0.09 ± 0.16 §, -0.05 < 
dns/dlnk < 0.02 §, £l < 0.05 and 12.5ei < r < 1 §. 
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IV. DISCUSSION 

The analysis of the general solution (|l^) confirms that 
the power-law inflationary scenario is just one among 
many suitable final stages for the inflaton dynamics as it 
was previously discussed in Ref. Moreover, as the ex- 
amples in the previous section indicate, a tensor to scalar 
ratio r converging to a constant value while the early uni- 
verse inflates is neither a necessary nor a sufficient con- 
dition for the corresponding inflationary model to yield 
perturbation spectra consistent with current CMB and 
LSS data in the range of observable scales. 

From these examples it can also be concluded that rad- 
ically different inflationary dynamics can yield the same 
asymptotic functional form for the tensor to scalar ratio. 
Therefore, the derivation of the inflaton potential from 
the information on the functional form of r will hardly 
be unique. Here it is important to note that by 'different 
dynamics' it is implied solution behaviors, which can be 
quite different for several values of the involved param- 
eters, even if the corresponding analytical expression for 
the potential, with unvalued parameters, is unique. 

The above conclusion must be regarded together with 
the fact that it seems not possible to obtain information 
on the inflaton potential beyond the exponential form, 
using exclusively the observational information on values 
of the spectral indices (ns and nr) and the running of 
the scalar index (dns / 'dlnfc) evaluated at the pivot scale, 
or on a constant central value of the tensor to scalar ratio 
. This would be a serious handicap for any program of 
reconstruction of the inflaton potential. 



A way to improve the situation described in the above 
paragraphs, could be to combine the information on An 
(the difference between the spectral indices) and the 
value of r, with the two first horizon-flow functions. 
An can be derived from observations and, according to 
Eq. (||), gives information on the scale dependence of r. 
With this information evaluated at the pivot scale, the 
underlying inflaton potential can be promptly sketched 
substituting expression ( p8| ) into Eq. (15). Moreover, 
here all three cases, power-law inflation, B ^ (if imitat- 
ing power-law inflation) and B = can be differentiated. 
In the case of power-law inflation, An and £2 will be zero. 
If B 7^ then, as discussed in the previous section, its 
absolute value must be extraordinarily small what, after 
ten or more e-foldings, leads to £3 — 1/C » 1.3712, a 
distinctively large, hardly suitable, value which will be 
the dominant contribution to An according to Eq. (|^). 
Finally, for B = 0, e 2 = CrAn = A, e 3 = and the 
potential will be given by Eq. (^9|). As shown in Fig. [l], 
realizations of this potential resemble the cases of mono- 
mial potentials with even order (Vt, £2 < 0), and inflation 
near a maximum (V77, £2 > 0) (see Ref. for exam- 
ples of such inflationary scenarios) allowing, therefore, 
to observe features of the inflaton potential beyond the 
exponential form characteristic of power-law inflation. 
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